Tensor-Structured Factorized Calculation of Two-Electron Integrals in a General Basis

In this paper, the problem of efficient grid-based computation of the two-electron integrals (TEI) in a general basis is considered. We introduce the novel multiple tensor factorizations of the TEI unfolding matrix which decrease the computational demands for the evaluation of TEI in several aspects. Using the reduced higher-order SVD the redundancy-free product-basis set is constructed that diminishes dramatically the initial number O(N2 b ) of 3D convolutions, defined over cross-products of Nb basis functions, to O(Nb) scaling. The tensor-structured numerical integration with the 3D Newton convolving kernel is performed in 1D complexity, thus enabling high resolution over fine 3D Cartesian grids. Furthermore, using the quantized approximation of long vectors ensures the logarithmic storage complexity in the grid-size. Finally, we present and analyze two approaches to compute the Cholesky decomposition of TEI matrix based on two types of precomputed factorizations. We show that further compression is possible via columnwise quantization of the Cholesky factors. Our “black-box” approach essentially relaxes limitations on the traditional Gaussian-type basis sets, giving an alternative choice of rather general low-rank basis functions represented only by their 1D samplings on a tensor grid. Numerical tests for some moderate size compact molecules demonstrate the expected asymptotic performance.